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Abstract 

The periodic Lorentz gas is the dynamical system corresponding to the free motion of a point particle in a periodic 
system of fixed spherical obstacles of radius r centered at the integer points, assuming all collisions of the particle 
with the obstacles to be elastic. In this Note, we study this motion on time intervals of order 1/r as r ^ 0^. 

Resume 

La limite de Boltzmann-Grad du gaz de Lorentz periodique en dimension deux d'espace. Le gaz 

de Lorentz periodique est le systeme dynamique correspondant au mouvement libre dans le plan d'une parti- 
cule ponctuelle rebondissant de maniere elastique sur un systeme de disques de rayon r centres aux points de 
coordonnees entieres. On etudie ce mouvement pour r ^ 0"*" sur des temps de I'ordre de 1/r. 



Version frangaise abregee 

On appelle gaz de Lorentz le systeme dynamique correspondant au mouvement libre d'une particule 
ponctuelle dans un systeme d'obstacles circulaires de rayon r centres aux sommets d'un reseau de R^, 
supposant que les collisions entre la particule et les obstacles sont parfaitement elastiques. Les trajectoires 
de la particule sont alors donnees par les formules (2). La limite de Boltzmann-Grad pour le gaz de Lorentz 
consiste a supposer que le rayon des obstacles r ^ 0"*", et a observer la dynamique de la particule sur 
des plages de temps longues, de I'ordre de 1/r — voir (3) pour la loi d'echelle de Boltzmann-Grad en 
dimension 2. 

Or les trajectoires de la particule s'expriment en fonction de I'application de transfert d'obstacle a 
obstacle Tr definie par (8) — oii la notation Y designe la transformation inverse de (7) — application 
qui associe, a tout parametre d'impact /i' G [—1,1] correspondant a une particule quittant la surface d'un 
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obstacle dans la direction w € S , le parametre d'impact /i a la collision suivante, ainsi que le temps s 
s'ecoulant jusqu'a cette collision. (Pour une definition de la notion de parametre d'impact, voir (6).) 

On se ramene done a etudier le comportement limite de I'application de transfert Tr pour r ^ 0+. 
Proposition 0.1 Lorsque < uj2 < oji et a = ^ ^ Q, Vapplication de transfert est approchee d 
0{r^) pres par I'application Ta,b,q,n definie d la formule (14). Pour u) G quelconque, on se ramene 
au cas ci-dessus par la symetrie (15). 

Les parametres A, B, Q,N mod. 2 intervenant dans I'application de transfert asymptotique sont definis 
a partir du dcvcloppcmcnt cn fraction continue (9) de a par les formules (11) ct (12). 

On voit sur ces formules que les parametres A, B,Q,N mod. 2 sont des fonctions tres fortement oscil- 
lantes des variables w et r. II est done naturel de chercher le comportement limite de I'application de 
transfert T,. dans une topologic faiblc vis a vis de la dcpendance en la direction uj. On montrc ainsi que, 
pour tout h' e [—1, 1], la famille d'applications ui i-^ Tr{h',ui) converge au sens des mesures de Young 
(voir par exemple [8] p. 146-154 pour une definition de cette notion) lorsque r — > 0+ vers une mesure de 
probabilitc P(s, h\h')dsdh indepcndante de lo : 

Theoreme 0.2 Pour tout $ S Cc(R+x] - 1, 1[) et tout h' e] - 1, 1[, la limite (16) a lieu dans L~(Si) 
faible-* lorsque r — > 0+, ou la mesure de probabilite P{s, h\h')dsdh est I'image de la probabilite ^ definie 
dans (17) par I'application {A,B,Q,N) TA,B,Q,N{h') de la formule (14). -De plus, cette densite de 
probabilite de transition P{s,h\h') verifie les proprietes (18). 

Le theoreme ci-dessus est le resultat principal de cette Note : il montre que, dans la limite de Boltzmann- 
Grad, le transfert d'obstacle a obstacle est decrit de maniere naturelle par une densite de probabilite de 
transition P{s,h\h'), oii s est le laps de temps entre deux collisions successives avec les obstacles (dans 
I'echelle de temps de la limite de Boltzmann-Grad), h le parametre d'impact lors de la collision future et 
h' celui correspondant a la collision passcc. 

Le fait que la probabilite de transition P{s, h\h') soit independante de la direction suggere I'hypothese 
d'independance (H) des quantites A, B,Q,N mod. 2 correspondant a des collisions successives. 
Theoreme 0.3 Sous I'hypothese (H), pour toute densite de probabilite £ Cc(R^ x S^), la fonction 
de distribution fr = fr(t,x,uj) de la theorie cinetique, definie par (3) converge dans L°°(R-)- x x S^) 
vers la limite (22) lorsque r 0+, oii F est la solution du probleme de Cauchy (21) pose dans I'espace 
des phases etendu {x,lj, s, ft,) G R^ x S"*- x R+x] — 1, 1[. 

Dans le cas d'obstacles aleatoires independants et poissonniens, Gallavotti a montre que la limite de 
Boltzmann-Grad du gaz de Lorentz obeit a I'equation cinetique de Lorentz (4). Le cas periodique est 
absolument different : en se basant sur des estimations (cf. [3] ct [7]) du temps de sortie du domaine 
Zr defini dans (1), on demontre que la limite de Boltzmann-Grad du gaz de Lorentz periodique ne pent 
pas etre decrite par I'equation de Lorentz (4) sur I'espace des phases R^ x classique de la theorie 
cinetique : voir [6]. Si I'hypothese (H) ci-dessous etait verifiee, le modele cinetique (22) dans I'espace des 
phases etendu fournirait done I'equation devant remplacer I'equation cinetique classique de Lorentz (4) 
dans le cas periodique. 



1. The Lorentz gas 

The Lorentz gas is the dynamical system corresponding to the free motion of a single point particle in 
a periodic system of fixed spherical obstacles, assuming that collisions between the particle and any of 
the obstacles are elastic. Henceforth, we assume that the space dimension is 2 and that the obstacles are 
disks of radius r centered at each point of Z^. Hence the domain left free for particle motion is 

Zr = {x £ \ dist(a;, Z^) > r} , where it is assumed that < r < 5. (1) 
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Assuming that the particle moves at speed 1, its trajectory starting from x € Zr with velocity a; € at 
time f = is f I— > (X^, Slr){t', x, w) e x given by 



Xr{t) = ^r{t) and tlr{t) = whenever Xr{t) G Zr, 

Xr{t + 0) = Xr{t - 0) axid nr{t + 0) = TZ[Xr{t)]nrit - 0) whcnevcr (t - 0) € , 

denoting ' = ^ and 7l[Xr{t)] the specular reflection on dZr at the point Xr{t) = Xr{t ± 0). Assume that 
the initial position x and direction u of the particle are distributed in x with some probability 
density = P'^{x,uj), and define 

fr{t,x,ui) := f''"{rXr{—t/r;x,co),Clr{—t/r;x,ui)) whenever x G Zr- (3) 

Wc are concerned with the limit of /r as r ^ 0+ in some appropriate sense to be explained below. In the 
2-dimensional setting considered here, this is precisely the Boltzmann-Grad limit. 

In the case of a random (Poisson), instead of periodic, configuration of obstacles, Gallavotti [5] proved 
that the expectation of fr converges to the solution of the Lorentz kinetic equation for {x,oj) e x S^: 

{dt+LO- V.)f{t, x,u)= [ {fit, x,iO- 2{uj ■ n)n) - f{t, x, uj)){uj ■ n)+dn , f _ = . (4) 

SI 

In the case of a periodic distribution of obstacles, the Boltzmann-Grad limit of the Lorentz gas cannot 
be described by a transport equation as above: sec [6] for a complete proof, based on estimates on the 
free path length to be found in [3] and [7]. This limit involves instead a linear Boltzmann equation on 
an extended phase space with two new variables taking into account correlations between consecutive 
collisions with the obstacles that are an effect of periodicity: see Theorem 4.1. 



2. The transfer map 

Denote by the inward unit normal to Zr at the point x G dZr, consider 

Tf = {{x,L0)edZr XS^\ ±10-71^ >0}, (5) 

and let F^/Z^ be the quotient of under the action of by translation on the x variable. For 
{x,aj) G F+, let Tr{x,Lu) be the exit time from x in the direction lo and hr{x,uj) be the impact parameter: 

Tr{x,cu) = mi{t > 0\ X + tiv £ dZr} , and hr{x,cu) = sm{cj,nx) ■ (6) 

Obviously, the map 

r+/z2 9 {x,uj) ^ {hr{x,(j),(j) e] - l,l[xS^ (7) 

coordinatizes F+/Z^, and we henceforth denote Yr its inverse. 

For each r e]0, consider now the transfer map Tj. :] — 1, l[xS^ R!i_x] — 1, 1[ defined by 

Tr{h',oj) = (rT^(y^(/l^a;)),/l^(X^(Tr(y^(/^^w));y^(/l',a;)),0^(T^(F^(/l^a;));Fr(/l^w)))). (8) 

For a particle leaving the surface of an obstacle in the direction uj with impact parameter h', the 
transition map Tr{h',u}) = {s,h) gives the (rescaled) distance s to the next collision, and the corresponding 
impact parameter h. Obviously, each trajectory (2) of the particle can be expressed in terms of the transfer 
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Figure 1. Left: the transfer map {s,h) = Tr{h',v), with h' = s'm{nj.i,v} and h = s'm{nx,v). Right: Particles leaving the 
surface of one obstacle will next collide with one of generically three obstacles. The figure explains the geometrical meaning 

of A,B,Q. 

map Tr and iterates thereof. The Boltzmann-Grad hmit of the periodic Lorcntz gas is therefore reduced 
to computing the hmiting behavior of Tr as r ^ and this is our main purpose in this Note. 

We first need some pieces of notation. Assume u = (wi, ^2) with < u)2 < wi, and a = ^2/^1 €]0, 1[\Q. 
Consider the continued fraction expansion of a: 



a= [0;ao,ai,a2,...] 



1 



(9) 



ao 



ai 



Define the sequences of convergents {pn,Qn)n>o and errors {dn)n>o by the recursion formulas 
Pn+i = anPn + Pn-i , = 1 , Pi = , rf„ = (-l)""^(g'„a - p„) , 

9n+l = anQn + Qn-l ?0 = , = 1 , 

and let 



N{a, r) = inf {n > 1 d„ < 2r-\/l + a2} , and k{a, r) 



2rv'l + - dM{a,r)-i 



dN(oL,r) 



(10) 



(11) 



Proposition 2.1 For each to = {cos9, sm9) with < 6 < j, set a = tan 6* and e = 2r^l + o? , and 

A{pi, r) = 1 ^ , B(a, r) = 1 , Q{a, r) = egjv(c«,r) • (12) 

In the limit r 0+, the transition map Tr defined in (8) is explicit in terms of A,B,Q,N up to 0{r^): 



Tr{h',Uj) = TA(a,r),B{a,r),Q(a,r),N{a,r){h') + (O(r^), 0) for each k' G] - 1, 1[. 

In the formula above 

TA,B,QMh') = {Q, h' - 2(-l)^(l - A)) if i-lfh' e]l - 2A, 1] , 

TA,B,QMh') = {Q', h' + 2(-l)^(l - B)) if (-1)^/1' e [-1, -1 + 2B[ , 

TA,B,Q,N{h') = {Q' + Q,h' + 2(-l)^(A - B)) if (-1)^/1' e [-1 + 2B,1- 2A] , 

for each {A,B,Q,N) e K :=]0, lpxZ/2Z, with the notation Q' = 
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(13) 



(14) 



y _ i-Q(i-B) 
i-A 



The proof uses the 3-term partition of the 2-torus defined in section 2 of [4] , following the work of [1] . 

For Lo = (cos0, sin^^) with arbitrary 6 G H, the map h' i— > Tr(h',uj) is computed using Proposition 2.1 
in the following manner. Set 9 — 9 ~ with m = [^{9 + |-)] and let uj = (cos 6', sin 6*). Then 

Tr{h' ,Lo) = {s,h) , where (s, sign(tan^)/i) = Tr(sign(tan^)/i', . (15) 



3. The Boltzmann-Grad limit of the transfer map Tj. 



The formulas (11) and (12) defining A, B, Q, N mod. 2 show that these quantities are strongly oscillating 
functions of the variables ui and r. In view of Proposition 2.1, one therefore expects the transfer map 
to have a limit as r ^ 0+ only in the weakest imaginable sense, i.e. in the sense of Young measures — 
see [8], pp. 146-154 for a definition of this notion of convergence. 

The main result in the present Note is the theorem below. It says that, for each h' G [—1,1], the family 
of maps Lu Tr{h',uj) converges as r ^ 0+ and in the sense of Young measures to some probability 
measure P{s, h\h')dsdh that is moreover independent of uj. 
Theorem 3.1 For each $ e Cc(R+ x [-1, 1]) and each h' 



CO 1 

^{Tr{h',-))-^ j J <P{s,h)P{s,h\h')dsdh in L°°{Sl) weak-* as r ^ 0+ , 



(16) 



-1 

where the transition probability P{s, h\h')dsdh is the image of the probability measure on K given by 

dii{A,B,Q,N) = -^lo<A<ilo<i;<i-Alo<Q<^--i--g ^^^^'^ {6n=o + Sn=i) (17) 
under the map K 3 [A, B, Q, N) i— > TA,B,Q,N{h') S R+ x [—1, 1]. Moreover, P satisfies: 

(s, h, h') (1 + s)P{s, h\h') is piecewise continuous and bounded on R+ x [—1, 1] x [—1, 1], 

(18) 

and P{s, h\h') = P{s, -h\ - h') for each h, h' G [-1, 1] and s > 0. 

The proof of (16-17) is based on the explicit representation of the transition map in Proposition 2.1 
together with Kloosterman sums techniques as in [2]. The explicit formula for the transition probability 
P is very complicated and we do not report it here, however it clearly entails the properties (18). 



4. The Boltzmann-Grad limit of the Lorentz gas dynamics 

For each r g]0, denote rf7+(a;, w) the probability measure on r+ that is proportional to w ■ Hxdxdu). 
This probability measure is invariant under the billiard map 

Br : r+ 9 {x,u}) I— > Br(a;, w) = (a; + Tr{x,u})ijO,TZ[x + Tr{x,u))oii]u}) G r+ . (19) 

For {x°,u°) e r+, set (a;",a;") = Bp(a;°,a;°) and a" = min(|a;f/a;J|, |a;JM|) for each n > 0, and define 

b'f = {A{an,r),B{an,r),Q{an,r),N{an,r) mod. 2) e K for each n > 0. (20) 
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We make the following asymptotic independence hypothesis: for each n > 1 and each ^ e C([— 1, 1] x K") 

1 

{H) ^lim y ^{hr,Uo,bl,...,b:;^)dj+{xo,Uo) = j ^ j ^ j ^{h',iVo,Pl,...,Pn)dfl{Pl)...dll{(}n) 



Under this assumption, the Boltzmann-Grad limit of the Lorentz gas is described by a kinetic model 
on the extended phase space R X X R+ X [—1,1] — unlike the Lorentz kinetic equation (4), that is 
set on the usual phase space R^ x S-*-. 

Theorem 4.1 Assume (H), and let /™ be any continuous, compactly supported probability density on 
R^ X . Denoting by R[6] the rotation of an angle 6, let F = F{t, x, w, s, h) be the solution of 



{dt+uJ-Vx-ds)F{t,x,u:,s,h) = j P{s,h\h')F{t,x,R[K - 2axcsm.{h')]uj,Q,h')dh' 

-1 

oo 1 

F{Q,x,w,s,h) =/'"(a;,w) j j P{T,h\h')dh'dT 



(21) 



where (a;, a;, s, h) runs through x S-^ x R;|J_x] — 1, 1[. Then the family (/r)o<r<^ defined in (3) satisfies 



oo 1 



fr-^ j j F{-,-, ■, s, h)dsdh in L°°(R+ x R^ x S^) weak-* as r ^ 0+. (22) 

-1 

For each (sq, ft-o) G 1^+ x [^1. 1], let (s„, hn)n>i be the Markov chain defined by the induction formula 

{Sni hn) = T/3^(/i„_i) for each n > 1 , together with a;„ = ^[2 arcsin(/i„_i) — 7r]a;„_i , (23) 

where /?„ G K arc independent random variables distributed under /i. The proof of Theorem 4.1 relies 
upon approximating the particle trajectory {Xr,^r){i) starting from {xq,uiq) in terms of the following 
jump process with values in R^ x S-^ x R_|_ x [—1, 1] with the help of Proposition 2.1 

(Xt,rit, S't,i/i)(xo,t^o, So, /lo) = (.-z^o + it^o,t^o, So - i,/io) for < t < sq , ^ ^ 

(24) 

{Xt,VLt, St, Ht){xQ,UjQ, 3q, ho) = {Xr„ + {t- S„)W„, W„, Sn+1 - t,hn) for Sn<t < S„+i . 

Unlike in the case of a random (Poisson) distribution of obstacles, the successive impact parameters on 
each particle path are not independent and uniformly distributed in the periodic case — likewise, the 
successive free path lengths on each particle path arc not independent with exponential distribution. The 
Markov chain (23) is introduced to handle precisely this difficulty. 
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